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ABSTRACT 


This  study  provides  an  aid  for  municipal  decision  makers 
of  the  Monterey  Peninsula.  Under  certain  assumptions,  the 
result  of  this  model  provides  optimal  annual  consumption 
given  the  current  water  supplies.  The  model  employs  present 
value  of  utility  of  private  consumption  by  peninsula  resi- 
dents as  the  objective  function.  Water  supplies  are  provided 
by  the  conjunctive  use  of  surface  reservoirs  and  a confined 
coastal  aquifer.  The  Carmel  Valley  watershed  was  modelled 
using  a statistical  application  of  the  Hydrologic  Balance 
equations  for  ground  water  systems. 
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I.  ORIGIN  OF  THE  PROBLEM  AND  THE  PROPOSED  MODEL 


The  following  study  was  motivated  by  the  1975-77  Cali- 
fornia drought  which  prompted  the  controversial  decision  in 
the  winter  of  1977  to  ration  the  consumption  of  water  to 
the  residents  of  the  Monterey  Peninsula.  Individual  consump- 
tion was  restricted  to  a monthly  maximum  not  to  exceed  a per 
capita  daily  consumption  rate  of  50  gallons.  By  comparison, 
average  annual  consumption  figures  for  the  normal  years 
immediately  preceeding  the  drought  reflected  a daily  rate  of 
135-150  gallons  per  person.1 

The  winter  of  1977-1978  brought  near  record  rainfalls  to 
the  Monterey  coast  that  ended  both  the  drought  and  the 
rationing  program  about  one  year  after  the  restrictions  were 
placed  in  effect. 

In  retrospect,  the  rationing  imposed  represented  a drastic 
response  to  the  water  shortage  (the  restriction  was  less 
than  l/3rd  of  normal  consumption  for  the  summer  months) . 

When  the  rains  came  in  December  of  1977,  existing  water 
supplies,  although  considerably  below  normal  levels  for  that 
time  of  year,  were  not  yet  considered  to  be  near  the  point 
of  exhaustion.  These  facts  considered,  a number  of  questions 


Unless  referenced  otherwise,  all  statistics  concerning 
the  Monterey  Peninsula  Water  Supply  and  consumption  were 
obtained  from  the  California  American  Water  Company. 
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arose  concerning  the  decisions  made.  First,  could  the  ration- 
ing have  been  made  less  drastic  by  enacting  milder  restric- 
tions? In  view  of  the  eventual  return  to  above  normal  rainfall 
conditions,  was  the  decision  made  too  hastily?  Finally,  if 
the  answer  to  either  of  the  preceeding  questions  is  yea,  then 
does  there  exist  some  procedure  that  the  decision-maker  can 
use  that  will  provide  him  with  a realistic  rationing  policy 
for  water  consumption  that  is  both  timely  and  optimal  with 
respect  to  some  reasonable  criterion. 

An  answer  to  the  last  question  is  the  subject  of  this 
study.  The  model  presented,  although  designed  specifically 
to  accomodate  the  Monterey  Peninsula  situation,  can  and  has 
been  generalized  to  fit  many  other  water  resource  allocation 
schemes,  particularly  in  the  area  of  flood  control  and/or 
multiple  use  water  resource  systems. 

The  approach  to  the  problem  was  suggested  by  Dr.  Alan 
Washburn  at  NPS,  Monterey.  The  stochastic  nature  of  rainfall, 
the  sole  replenishment  source  of  water  for  the  peninsula, 
suggests  some  form  of  a Markov  decision  process  since  subse- 
quent supplies  of  water  are  affected  by  both  deterministic 
(consumption)  and  random  (rainfall)  processes. 

Utility  of  water  to  the  population  was  selected  as  the 

measure  of  benefit  to  be  used  in  the  objective  function 

2 

because  of  the  life  style  of  the  Peninsula.  Approximately 


2 

The  term  "Monterey  Peninsula"  shall  be  used  to  represent  an 
area  known  to  the  California  State  Department  of  Water  Resources  as 
Zone  11  and  includes  Carmel  Valley,  Carmel,  Monterey  and  parts 
of  Seaside. 


86%  of  the  water  metered  to  the  peninsula  is  used  to  meet 
residential  and  other  human  uses.  Less  than  1.3%  is  used 
in  local  industry.  For  this  reason  it  is  difficult  to  attach 
economic  values  to  water  that  other  models  of  this  type  have 
done.  Whether  the  actual  utility  of  water  for  the  population 
can  be  measured  is  a question  that  will  be  discussed  later 
in  detail. 

Assuming  that  such  a utility  function,  U(D),  can  be 
determined  for  the  population  where  D = water  available  for 
consumption,  then  a natural  objective  of  the  decision-maker 
is  to  establish  a consumption  policy  that  will  maximize  the 
total  utility  of  consumption. 

The  nature  of  the  hydrological  cycle  of  the  Carmel  Valley 
watershed3  relies  on  yearly  rainfall  to  provide  the  runoff 
necessary  to  replenish  water  supplies  depleted  by  consumption 
and  other  losses.  Therefore,  the  supplies  available  for  future 
consumption  are  dependent  upon  both  current  consumption  and 
future  runoff.  The  current  or  near-future  consumption  is 
determined  by  policy  while  future  runoff  is  dependent  on 
rainfall  which  is  random. 

It  is  reasonable  to  assume  that  the  population  derives 
utility  from  future  as  well  as  present  consumption.  It  is 
equally  reasonable  to  assume  that  the  utility  of  future 


3The  255  square  miles  of  the  upper  and  lower  Carmel 
Valley  constitutes  the  major  source  of  rainfall  runoff  used 
to  recharge  the  Monterey  Peninsula  water  supply. 
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consumption  is  not  as  important  as  that  of  present  or  near 
future  consumption.  For  this  reason,  the  concept  of  present 
value  of  the  utility  of  future  supplies  is  useful  and  the 
model  in  fact  employs  just  such  a concept  in  its  objective 
function. 

The  objective  then  is  to  select  a policy  for  current  or 
near  future  consumption  that  will  maximize  the  total  utility 
associated  with  that  consumption  plus  the  present  value  of 
the  utility  of  all  future  consumption,  the  supply  for  which 
is  not  yet  known  by  the  policy  maker.  What  is  known  once 
consumption  policy  is  set  is  the  amount  of  water  lost  from 
storage.  Now,  if  rainfall  could  be  predicted  and  its  effect 
on  water  supplies  known,  the  objective  function  would  be  com- 
plete. However,  for  this  model,  future  rainfall  cannot  be 
predicted.  What  can  be  estimated  is  the  probability  distri- 
bution of  annual  rainfall  which  can  then  be  used  to  compute 
expected  utility  of  future  consumption. 

The  final  objective  of  the  model,  therefore,  is  to  choose 
a consumption  decision  from  among  all  possible  decisions 
available  that  will  maximize  the  sum  of  the  utility  of  that 
decision  plus  the  present  value  of  the  expected  value  of  the 
utility  of  future  consumption. 

The  decision  period  selected  is  one  year  in  length.  Since 
timeliness  of  the  decision  is  a factor,  the  date  selected  for 
the  model's  decision  is  the  first  of  May.  This  is  a particu- 
larly meaningful  date  since  it  represents  not  only  the  end  of 
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annual  rainfall  for  seven  months  but  also  the  beginning  of 

the  period  of  heaviest  daily  consumption  of  water. 

Graphically  the  model  can  be  represented  as  a series  of 

stages  (see  figure  1) , one  for  each  year  into  the  future. 

Entering  the  nth  stage  are  the  state,  Xn>  representing  the 

amount  of  water  stored  in  the  supply  system,  a decision  Dn 

representing  allowable  consumption  and  an  amount  of  runoff, 

R . Leaving  the  stage  is  the  new  state  X , that  will  enter 
n n- 1 

the  (n-l)th  stage.  Finally,  associated  with  each  stage  is 
a value  U(Dn)  that  represents  the  utility  associated  with  Dn- 


, 4 

The  concept  of  the  model  is  now  complete.  The  following 
are  required  for  the  model  to  work: 


A 

The  final  model  required  is  somewhat  more  complex  as  shuwn 
in  Section  III.  A complete  mathematical  description  of  the 
model  is  given  in  Appendix  A,  containing  final  parameter  and 
function  descriptions. 

■ 
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1.  The  utility  function  U (D)  which  associates  an  arbitrary 
value  U(D)  with  the  consumption  quantity  of  water,  D. 

2.  The  replenishment  function  that  will  estimate,  for 
a given  rainfall,  the  amount  of  water  available  for  return 
to  the  water  supply. 

3.  The  maximum  total  storage  capacity  of  the  system. 

4.  The  probability  distribution  for  rainfall. 

Each  of  these  sub-models  is  developed  in  one  of  the 

remaining  sections  of  this  study  followed  by  a discussion  of 
the  results  of  the  model  under  various  assumptions.  The  study 
concludes  with  some  additional  remarks  and  recommendations 
for  future  study  that  resulted  from  the  research  associated 
with  this  study. 

The  majority  of  hydrological  information  and  data  was 
provided  by  the  Monterey  County  Flood  Control  and  Water 
Conservation  District  and  the  California  American  Water 
Company.  Their  eager  assistance  and  full  cooperation  have 
been  invaluable  in  the  completion  of  this  project. 


II.  RAINFALL  - THE  SOLE  SOURCE  OF  WATER 


The  most  extensive  rainfall  records  have  been  maintained 
at  the  Forest  Lake  holding  reservoir  near  the  coast.  Records 
there  have  been  kept  since  1896  with  only  a few  missing 
values.  Additional  rainfall  data  was  examined  at  three  other 
sites;  Pacific  Grove,  Los  Padres  Reservoir  and  San  Clemente 
Reservoir.  The  Pacific  Grove  and  Los  Padres  records  have 
been  kept  since  the  mid  1940' s while  San  Clemente  records  reach 
back  to  1921. 

The  rainfall  in  the  upper  Carmel  Valley  appears  to  be 
about  5 inches  more  annually  than  on  the  coast.  This  fact 
is  consistent  with  the  climatological  conditions  of  the  cen- 
tral California  coast.  Of  the  four  sites,  the  Los  Padres 
dam  rainfall  data  seems  to  have  the  highest  correlation  to 
the  runoff  data  computed  in  Section  III. 

Reference  1 presents  some  analytical  procedures  for 
accurately  modelling  the  total  precipitation  volume  for  a 
watershed  using  the  rainfall  data  from  several  gauge  sites. 

One  method  employs  polygonal  areas  of  constant  precipitation. 
Another  uses  contour  lines  of  equal  precipitation.  All  of 
them,  however,  seem  to  rely  on  a reasonably  dense  distribu- 
tion of  rain  gauge  data  across  the  watershed.  Such  data  is 
not  presently  available  for  the  Carmel  Valley. 

The  data  that  is  available  suggests  that  use  of  single 
site  rainfall  data  could  provide  an  adequate  indicator  of 


total  precipitation  over  the  entire  watershed.  The  simple 
correlations  between  the  four  sites  examined  indicate  a 
strong  dependence.  The  small  areal  extent  as  well  as  the 
local  climate  seem  to  preclude  independent  rainfall  behavior 
within  the  watershed. 


Since  Los  Padres  rainfall  provided  the  best  model  of 
runoff,  it  was  selected  as  the  annual  rainfall  to  be  used 
in  computation  of  expected  utilities.  Successive  annual 
rainfall  measurements  at  Los  Padres  were  assumed  to  be 
independent  random  variables.  Serial  correlation  of  existing 
data  was  low  and  its  possible  effects  were  ignored. 

A cumulative  probability  distribution  (CDF)  of  annual 
rainfall  was  then  constructed  from  the  existing  30  data  points 
and  is  shown  in  Figure  2.  The  CDF  can  now  be  used  to  graphically 
extract  the  probability  mass  function  required  by  the  model. 

Table  I contains  the  probability  mass  function  used  by  the 
final  model. 


PADRES  RAINFALL  (INCHES) 


Rainfall 

(inches) 

Less  than  1. 

1.5- 2. 5 

2. 5- 3. 5 

3. 5- 4. 5 

4. 5- 5. 5 

5. 5- 6. 5 

6. 5- 7. 5 

7. 5- 8. 5 

8. 5- 9. 5 

9.5- 10.5 

10.5- 11.5 

11. 5- 12.5 

12.5- 13.5 

13.5- 14.5 

14.5- 15.5 

15.5- 16.5 

16.5- 17.5 

17.5- 18.5 


Probability 
of  Occurrence 

5 .002 

.0005 
.0005 
.0005 
.005 
.006 
.006 
.006 
.009 
.011 
.011 
.015 
.016 
.024 
.035 
.047 
.078 
.102 


Rainfall 

(inches) 

18.5- 19.5 

19.5- 20.5 

20.5- 21.5 

21.5- 22.5 

22.5- 23.5 

23.5- 24.5 

24.5- 25.5 

25.5- 26.5 

26.5- 27.5 

27.5- 28.5 

28.5- 29.5 

29.5- 30.5 

30.5- 31.5 

31.5- 32.5 

32.5- 33.5 

33.5- 34.5 

34.5- 35.5 
Greater  than 


Probability 
of  Occurrence 

.12 

.06 

• 0,4 

.035 

.025 

.02 

.02 

.015 

.015 

.015 

.015 

.015 

.013 

.012 

.011 

.010 

.010 

.5  .1745 


Table  I 


III.  THE  MONTEREY  PENINSULA  WATERSHED  - 
A STUDY  IN  HYDROLOGY 


A major  effort  of  this  study  involved  the  question  of 
replenishment  of  water  supplies  to  the  storage  system  given 
the  current  years  rainfall.  This  is  actually  a generaliza- 
tion of  several  more  specific  questions. 

How  much  storage  is  available  for  collecting  the  water 
generated  by  rainfall?  How  much  is  lost  to  the  system  and 
how?  Does  runoff  due  to  rainfall  depend  only  upon  current 
rainfall  or  are  there  other  factors  to  consider?  The  answers 
to  these  questions  provide  a model  of  the  watershed  that  is 
sufficient  to  complete  the  general  model. 

A.  TOTAL  CAPACITY  OF  THE  SYSTEM 

Recent  studies  [Ref.  2, 3, 4, 5]  have  been  made  of  the  Carmel 
Valley  watershed,  the  major  source  of  surface  and  ground 
supplies  of  usable  water  for  the  peninsula.^  One  of  these 
studies  was  conducted  in  1974  by  the  California  Department 
of  Water  Resources  with  a revision  to  the  study  submitted  in 
December  1977.  Its  purpose  was  to  "conduct  an  independent 
study  of  the  yield  of  the  ground  water  basins  in  the  service 
area  ..."  (of  Zone  11)  [Ref.  2].  Zone  11  comprises  the  entire 


^There  are  wells  in  the  Seaside  area  that  provide  limited 
supplies  but  are  not  considered  part  of  the  Carmel  Valley 
watershed.  Their  input  to  the  system  is  small  and  will  be 
treated  as  constant. 


Carmel  Valley,  Monterey  Peninsula  and  Seaside  area.  The 
study  addressed  surface  storage  and  two  principal  aquifers, 
the  Carmel  River  Valley  alluvium  and  the  Seaside  aquifer. 


r 
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; 


\ 
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The  results  of  the  study  have  provided  much  of  the  infor- 
mation needed  for  establishing  the  storage  capacity  of  Jt he 
entire  system. 

The  Carmel  Valley  watershed  drains  about  255  square  miles 
of  terrain.  The  boundaries  are  considered  impervious  to 
flow,  thus  creating  an  enclosed  basin  whose  only  subsurface 
outlet  for  flow  occurs  at  a narrow  opening  or  crack  containing 
porous  alluvium  below  the  mouth  of  the  Carmel  River. 

There  are  two  man-made  reservoirs  located  at  the  Los  Padres 
and  San  Clemente  dams.  The  surface  water  capacity  of  the  Los 
Padres  reservoir  is  estimated  at  3000  acre  feet  with  a maximum 
surface  area  of  67  acres.  The  San  Clemente  reservoir  surface 
storage  is  estimated  at  2154  acre  feet  with  a maximum  surface 
area  of  53  acres  [Ref.  2]. 

The  Los  Padres  dam  is  located  approximately  6 miles  up- 
stream from  the  San  Clemente  dam.  Levels  at  San  Clemente 
reservoir  are  controlled  by  the  regulation  of  spillage  from 
Los  Padres  dam  and  by  diversions  from  the  reservoir  by  the 
California  American  Water  Company  to  meet  consumption  demands. 

The  drainage  area  above  the  San  Clemente  dam  is  about  125 
square  miles  and  although  the  upper  subsurface  is  porous, 
there  are  no  ground  water  basins.  The  lower  Carmel  Valley 
contains  the  remaining  130  square  miles  of  drainage  area,  the 
runoff  from  which  drains  primarily  into  the  Carmel  River. 
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By  far  the  major  storage  capacity  for  consumptive  water 


is  the  Carmel  Valley  alluvium.  It  is  estimated  to  extend 
13  miles  from  the  mouth  of  the  Carmel  River  to  a point  about 
1 mile  south  of  the  juncture  of  the  Tularcitos  creek  and  the 
Carmel  river.  According  to  the  Dames  and  Moore  study,  / 

The  alluvium  beneath  the  floor  of  the  Carmel 
Valley  ranges  in  width  at  its  surface  expression 
from  500  feet  in  upper  valley  areas  to  approxi- 
mately 1 mile  at  the  coast.  It  is  generally  less 
than  one  half  mile  wide.  Its  depth  ranges  from 
about  50  feet  in  upper  valley  areas  to  more  than 
150  feet  in  the  lower  valley  near  the  coast.  [Ref.  3] 

The  areal  boundaries  of  the  alluvium  are  estimated  to 

include  4210  acres  at  the  surface.  The  alluvium  consists 

primarily  of  "sand  and  gravel,  much  of  the  gravel  being  boulder" 

[Ref.  2].  The  results  of  the  1974  study  suggest  a total 

volume  of  the  alluvium  of  321,900  acre  feet  with  an  average 

thickness  of  76.5  feet. 

Individual  drillers  logs  at  various  well  sites  were 
analyzed  for  specific  yield6  of  the  materials  extracted. 

The  specific  yields  ranged  from  .1946  to  .2749  with  a mean 
value  of  .2359.  Thus  the  volume  of  extractable  water  during 
peak  storage  in  the  Carmel  Valley  aquifer  can  be  estimated  at 
321,900  x .2359  = 75,936  acre  feet.  In  December  1977  this 
estimate  was  revised  upward  for  the  western  portion  of  the 
aquifer  by  the  State  Water  Resources  Board  [Ref.  5]. 

6Specific  Yield  = % of  the  total  volume  of  the  alluvium 
occupied  by  the  ultimate  volume  of  water  that  can  be  released 
from  the  saturated  alluvium. 
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Although  the  report  acknowledges  over  78,000  acre  feet 
of  water  that  can  be  drawn  from  the  full  aquifer,  this  model 
will  be  based  on  the  safe  yield  concept  discussed  by  Hall  & 

Dracup  [Ref.  6].  Safe  yield  will  refer  to  "that  amount  of 
water  which  can  be  withdrawn  annually  without  causing  an 
undesired  influence  in  the  basin"  [Ref.  6].  There  are  three 
major  factors  that  need  to  be  considered  when  determining 
safe  yield. 

First,  when  an  aquifer  is  drawn  excessively  low,  the 
reduction  in  pressure  of  the  water  places  a greater  burden 
of  the  weight  of  the  alluvium  upon  the  alluvium  itself.  This 
causes  grains  and  pebbles  of  the  alluvium  to  shift  resulting 
in  tighter  structure,  reduced  specific  yield  and  less  storage 
capacity.  This  phenomenon  is  known  as  land  subsidence  and 
is  not  considered  a factor  in  the  alluvium  below  the  Carmel 
riverbed. 

Second,  long  term  use  of  an  aquifer  requires  annual 
recharge  from  percolation  of  runoff  due  to  rainfall.  If 
annual  extractions  exceed  the  long  term  annual  recharge 
capability  of  runoff,  then  over  time  the  water  table  will  be 
lowered.  This  phenomenom  has  caused  damage  to  many  valuable 
ground  water  tables  in  the  past  and  only  in  recent  years  has 
there  been  careful  monitoring  of  the  long  term  balance 
between  extraction  and  recharge. 

The  recharge  rate  of  the  Carmel  Valley  alluvium  is  remarka- 
bly rapid  due  to  its  high  permeability.  Recharge  is  accomplished 

1 
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by  percolation  of  stream  flow  of  the  Carmel  river  through  the 
riverbed.  Dames  and  Moore  [Ref.  3]  place  percolation  rates 
of  the  streambed  as  high  as  100  cubic  feet  per  square  feet 
per  second,  while  the  Zone  11  study  observed  a recovery  rate 
of  "8  to  23  feet  per  month  in  November  1972  following  the  last 
dry  period"  [Ref.  5]. 

The  50  year  mean  annual  full  natural  flow  of  the  Carmel 
river  at  the  San  Clemente  dam  is  estimated  to  be  61,900  acre 
feet.  Consequently,  the  safe  yield  for  maintenance  of  the 
water  table  is  quite  high,  perhaps  exceeding  50,000  acre  feet. 

Third,  coastal  aquifers  face  a unique  problem  when  the 
alluvium  extends  beyond  the  seashore.  In  effect  the  alluvium 
is  saturated  below  the  land  side  with  fresh  water  and  below 
the  ocean  side  with  salt  water.  Within  the  alluvium  there 
exists  a boundary  or  face  consisting  of  a mix  of  salt  and 
fresh  water.  This  boundary  is  maintained  close  to  the  shore 
line  by  the  pressure  of  the  water  table  which  is  elevated 
above  sea  level,  producing  a seaward  flow  of  fresh  water. 

Sea  water  is  heavier  than  fresh  water.  Should  the  height 
of  the  water  table  be  reduced  below  certain  levels,  the  heavier 
sea  water  could  set  up  reverse  flow  landward  beneath  the  fresh 
ground  water.  Since  commercial  wells  are  normally  drilled 
the  entire  depth  of  the  alluvium  for  maximum  withdrawal,  wells 
in  such  locations  would  no  longer  be  of  any  use  due  to  salt 
water  contamination  in  their  lower  parts  (see  Figure  3) . 
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There  is  a simple  thumb  rule  for  the  water  table  elevation 
required  to  prevent  sea  water  intrusion.  Consider  the  columns 
of  water  in  Figure  3 represented  by  H and  H+t.  For  H+t 
to  be  fresh  water,  it  must  be  equal  in  weight  to  the  column  H 
of  the  heavier  sea  water.  Using  specific  gravities  of  ,1  and 
1.025  for  fresh  and  sea  water,  the  following  relationship 
holds: 

(H+t)  = 1.025  H , 


or  H = 40t,  or  for  every  foot  that  the  water  table  lies  above 
sea  level,  40  feet  of  fresh  water  can  be  maintained  below 
sea  level. 

The  deepest  portion  of  the  Carmel  Valley  aquifer  lies 
next  to  the  coast.  Well  depth  measurements  indicate  a maxi- 
mum depth  of  160  feet  below  sea  level  to  impervious  strata. 

The  water  table  at  the  coastline,  therefore,  must  be  main- 
tained at  or  above  4 feet  above  sea  level  to  prevent  sea 
water  intrusion.  This  is  by  far  the  most  restrictive  require- 
ment under  the  safe  yield  concept  discussed  earlier.  That 
is,  this  model  will  consider  the  gross  capacity  of  the  aquifer 
to  be  that  yield  of  the  alluvium  that  lies  above  the  minimum 
water  table  elevation  required  to  preclude  sea  water  intrusion, 
under  present  methods  of  extraction.7 


7There  exist  methods  (most  quite  expensive)  which  enable 
inland  lowering  of  a coastal  water  table  while  maintaining 
a pressure  "barrier"  at  the  coast  line  to  prevent  sea  water 
intrusion. 
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This  approach  is  considered  quite  conservative  by  some 
hydrologists  as  illustrated  by  the  following  comment  by  an 
engineer  with  the  California  State  Water  Resources  board. 

Pumping  the  ground  water  until  its  levels  are 
below  sea  level  will  not  cause  sea  water  intru- 
sion as  long  as  there  exists  a mound  of  fresh 
water  with  an  elevation  greater  than  sea  level 
between  the  ocean  and  the  area  pumped.  This 
could  be  done  on  a short  time  basis  without 
injection  barriers  [Ref.  7]. 

Another  report  on  ground  water  resources  of  the  Carmel 
Valley  was  submitted  by  the  Carmel  Valley  Master  Plan  Study 
Committee  (Ref.  4].  It  used  the  data  from  the  Zone  11  study 
to  construct  cross  sectional  views  of  the  alluvium  in  an 
effort  to  accurately  determine  what  portion  of  the  alluvium 
lay  sufficiently  above  sea  level  to  preclude  sea  water  intru- 
sion. The  results  were  an  alluvial  volume  of  116,214  acre 
feet  which,  given  a permeability  of  .2359,  computes  to  27,500 
acre  feet  of  storage. 

Total  usable  storage  of  the  system,  therefore,  is  27,500 
acre  feet  of  ground  storage  plus  5,000  acre  feet  of  surface 
storage  in  the  reservoirs.  This  figure  should  be  considered 
the  most  conservative  estimate  of  maximum  storage  capacity 
of  the  watershed.  There  are  many  who  claim  the  figure  is 
considerably  higher  while  there  are  still  some  who  doubt  it 
is  that  high.  Some  of  the  differences  of  opinion  among  the 
"experts"  seem  to  be  linked  to  differing  points  of  view  on 
what  future  growth  should  be  permitted  on  the  Monterey 
Peninsula,  an  issue  totally  irrelevant  to  this  study. 


This  study  will  use  33,000  acre  feet  as  the  maximum 
static  capacity  for  usable  storage  of  water  by  the  Monterey 
Peninsula  watershed.  It  is  static  in  the  sense  that  this 
much  water  is  extractable  even  after  all  runoff  has  ceased. 

In  the  spring  of  the  year  when  this  model  considers 

alternative  decisions,  the  storage  system  may  not  be  static, 

particularly  if  the  previous  rains  were  at  or  above  normal 

levels.  It  is  possible  after  heavy  rainfall  for  stream  flow 

to  continue  into  the  reservoirs  and  the  Carmel  riverbed  as 

s t 

late  as  July.  If  the  static  storage  on  May  1 is  at  maximum 
capacity,  then  there  exists  a period  of  simultaneous  consump- 
tion and  replenishment  until  all  flows  cease.  This  additional 
source  of  residual  runoff  after  complete  recharge  will  be 
modelled  as  "virtual"  capacity  above  and  beyond  the  maximum 
static  capacity  already  discussed.  The  motivation  for 
modelling  it  in  this  manner  will  be  clarified  in  the  following 
section. 

B.  REPLENISHMENT  - A STATISTICAL  MODEL  OF  RUNOFF 

Several  references  texts  and  articles  [Ref.  1,8,9,10,11, 

12]  were  consulted  for  replenishment  modelling  techniques. 

Most  of  the  hydrological  models  discussed  utilized  physical 
models  of  runoff.  The  physical  characteristics  of  the  alluvium 
are  carefully  studied,  permeabilities  established,  runoff 
characteristics  of  the  terrain  are  estimated,  climatological 
data  collected  for  coefficients  of  evaporation,  types  and 
density  of  ground  vegetation  analyzed  to  determine 
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evapotranspiration  rates,  suburban  and  urban  development 
studied  for  their  associated  runoff  characteristics,  storm 
intensities  estimated  and  their  frequencies  recorded,  coeffi- 
cients of  storage,  permeability  and  transmissibility  estimated 
and  many  other  physical  properties  explored. 

From  such  data  parameters  are  established  that  enable  the 
hydrologist  to  simulate  the  behavior  of  the  watershed  using 
one  of  a variety  of  mathematical,  physical  and/or  analog 
type  models  depending  upon  his  objective.  The  parameters 
necessary  for  this  type  of  modelling  of  the  peninsula  water- 
shed are  not  available  nor  are  they  considered  obtainable 
within  the  scope  of  this  study. 

More  recently,  studies  in  hydrology  have  applied  statisti- 
cal approaches  to  runoff  models,  particularly  when  the  histori- 
cal data  base  is  available  and  the  physical  characteristics 
are  not  well  known.  The  technique  employs  a fundamental 
relationship  known  as  Hydrologic  Balance  which  takes  the 
following  form: 


INFLOW  - OUTFLOW  = CHANGE  IN  STORAGE 


Elements  of  inflow  and  outflow  are  illustrated  in  Figure 
4.  According  to  Peters  [Ref.  1],  the  specific  basin  balance 
may  contain  any  or  all  of  the  following: 
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INFLOW 


OUTFLOW 


CHANGE  IN  STORAGE 


Precipitation 

Surface  stream 
inflow 

Imported  water 


Imported  sewage 
Subsurface  inflow 
Desalinated  water 


Consumptive  Use 
Evaporation 


Exported  fresh 
water 


Reservoirs 
Ground  water 

Saline  water  intrusion 


Exported  sewage  Subsidence 

Subsurface  outflow  Perched  water 

Surface  stream 
outflow 

for  irrigation 
for  artificial 
recharge 


Agricultural  drains 


For  the  Monterey  Peninsula  watershed,  the  Hydrologic 
Balance  was  used  to  obtain  accurate  estimates  of  total  inflow. 
Since  all  inflow  is  due  ultimately  to  annual  rainfall,  total 
annual  inflow  data  could  then  be  compared  to  annual  rainfall 
data  in  hopes  of  determining  a function  for  estimating  the 
unknown  inflow  as  a result  of  known  rainfall. 

The  following  statistics  were  available  from  the  indicated 
sources : 

1.  Monthly  flows  of  the  Carmel  River  at  Carmel  from 
annual  reports  of  the  U.S.  Geological  Survey  [Ref.  13]. 

2.  Rate  of  subsurface  flow  into  the  ocean  as  esti- 
mated by  the  California  State  Water  Resources  Board  [Ref.  2]. 

3.  Total  diversions  from  the  Los  Padres/San  Clemente 
Reservoirs  from  the  records  of  the  California  American  Water 
Company  (Cal  Am) . 
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4.  Total  monthly  extractions  from  the  Carmel  Valley 
aquifer  by  Cal  Am  as  stated  in  their  records. 

5.  Annual  private  consumption  from  the  Carmel  Valley 
aquifer  as  estimated  by  the  California  State  Water  Resources 
Board. 

6.  Monthly  reservoir  levels  recorded  by  Cal  Am. 

7.  Annual  average  change  in  well  levels  in  the  Carmel 
Valley  aquifer  from  the  records  of  the  Monterey  Flood  Control 
and  Water  Conservation  District  [Ref.  14). 

The  data  was  collected  by  water  year  defined  as  1 October 
through  30  September  in  order  to  insure  as  static  a reading 
of  the  system  as  possible.  In  the  Carmel  Valley,  rainfall 
occurs  between  the  months  of  October  and  April  with  little 
or  no  rain  throughout  the  summer,  so  that  by  1 October  all 
runoff  can  be  considered  complete.  This  provides  the  most 
accurate  computation  of  runoff  following  a winter  rain. 
Unfortunately  complete  records  have  been  kept  only  since  1963. 
They  are  shown  in  Table  II. 

The  computation  of  runoff  was  done  as  follows:  Runoff 
(Inflow)  = Change  in  storage  plus  Outflow,  where  Change  in 
storage  = Change  in  Reservoirs  plus  change  in  ground  water 
and  Outflow  = (Carmel  River  flow  at  Carmel)  plus  (subsurface 
flow  at  Carmel)  plus  (Reservoir  diversion  by  Cal  Am)  plus 
(Well  production  by  Cal  Am)  plus  (Private  Well  production) . 

The  resulting  runoff  computations  are  tabulated  in  Table  III  by 
year  with  the  year's  corresponding  rainfall  total  at  Los  Padres 


dam. 


Stream  Flow  Subsurface  Producticx 

at  River  Mouth  Flew  at  Mouth  by  Cal  Am 
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Exploratory  data  analysis  was  performed  on  the  runoff 
and  rainfall  data,  the  results  of  which  were  quite  promising. 
Figure  5 shows  a plot  of  Runoff  versus  Rainfall.  The  plot 
has  definite  structure  with  an  apparent  linear  relationship 
except  for  the  1969  runoff  figure  that  was  recorded  during 
flood  conditions.  That  value  was  subsequently  discarded  for 
two  reasons.  First,  it  is  highly  likely  that  it  is  the  least 
accurate  value  in  view  of  the  flooding  that  occurred  when 
measurements  were  taken.  Second,  the  primary  concern  is  for 
accurate  estimates  of  normal  to  below  normal  runoff  since  near 
flood  conditions  will  easily  recharge  the  system. 

A least  squares  fit  of  the  remaining  runoff  versus  annual 
rainfall  produces  a good  model  in  a statistical  sense;  however 
there  are  certain  qualities  of  the  runoff  data  revealed  in 
Figure  5 that  imply  additional  dependency.  All  the  circled 
points  which  lie  above  the  fitted  line  have  something  in 
common.  They  all  represent  values  that  occurred  in  a year 
following  a previous  year  of  above  average  runoff.  On  the 
other  hand  the  circled  points  lying  below  the  line  represent 
values  which  occurred  following  a year  of  below  average  runoff. 
This  suggests  that  previous  runoff  has  an  impact  on  current 
runoff. 

A multiple  regression  of  runoff  as  a function  of  rainfall 
and  the  one  year  lag  of  runoff  produced  slightly  better 
statistics  and  far  better  parallels  of  the  trends  of  the  actual 
data.  Unfortunately  the  decision  model's  state  vector  now 
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has  two  dimensions  (Rainfall,  Previous  Runoff) , each  of  which 
has  considerable  range.  A simple  heuristic  approach  was 
tried  to  reduce  the  number  of  states  required  overall. 

Assume  that  annual  runoff  depends  not  only  upon  the  amount 
of  annual  rainfall  but  also  upon  the  wetness  of  the  soil 
prior  to  the  rainy  season.  Furthermore,  assume  that  condition 
of  the  soil  can  be  categorized  in  broad  terms  as  dry,  normal 
or  wet  (called  wetness) . Now  the  second  dimension  of  the 
state  vector,  namely  wetness,  need  only  range  over  three 
values. 

The  following  heuristic  gave  the  best  results.  Let  wet- 
ness, Wn,  be  1,  2,  or  3 corresponding  to  dry,  normal  and  wet 
soil  conditions.  The  value  of  is  found  from  the  previous 
years  runoff  as  follows: 
r 

1 If  the  previous  runoff  was  less  than 
11,601  acre  feet  (one  standard  deviation 
below  the  mean®) 

Wetness  = / 3 If  the  previous  runoff  was  greater  than 

148793  acre  feet  (one  standard  deviation 
above  the  mean) 

2 Otherwise 

' v 

The  resulting  model  which  is  presented  below  far  exceeded 
expectations  in  accuracy  and  estimating  power: 


O 

The  mean  and  standard  deviation  referred  to  are  the 
sample  mean  and  deviation  computed  from  all  the  runoff  data 
available. 


32 


I 


Runoff  » 5296  (Rain)  + 8433  (Wetness)  - 71,059  Acre  feet 


except  that  in  no  case  is  runoff  estimated  to  be  less  than 
zero. 

Figure  6 shows  a plot  of  actual  Runoff  vs  the  estimated 
Runoff  by  the  model. 


C.  VIRTUAL  CAPACITY 

The  above  model  owes  much  of  its  accuracy  and  simplicity 
to  the  selection  of  October  as  the  start  of  the  annual  period, 
enabling  an  accurate  estimate  of  runoff  with  only  two  pieces 
of  information  required;  annual  rainfall  for  the  year  and 
wetness . 

Since  streamflow  can  continue  for  some  time  after  heavy 
rainfall  has  ceased,  it  is  possible  for  the  reservoirs  and 
the  alluvium  to  continue  recharging  for  some  time  after  the 
first  of  May,  the  date  marking  the  beginning  of  the  decision 
period.  What  this  amounts  to  is  additional  storage  of  sub- 
surface flows  that  have  not  yet  been  removed  by  stream  flow 
to  the  ocean.  This  additional  storage,  called  "virtual 
storage",  is  considered  available  only  during  May,  June  and 
July.  Therefore  the  most  that  it  could  ever  be  is  the  maximum 
expected  demand  from  the  system  during  these  months.  Histor- 
ically this  amounts  to  about  5,000  acre  feet. 
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(1000  ACRE  FEET) 


D 


SUMMARY  OF  THE  MONTEREY  PENINSULA  WATERSHED  MODEL 


1.  Capacity 

1 May  Maximum  Capacity  = Static  Capacity  + Virtual  Capacity 

= 33,000  + 5,000 

= 38,000  acre  feet 


2.  Replenishment 

Given  R inches  of  rain  and  wetness  = W, 

Total  Annual  Runoff  = 5296*R  + 8433»W  - 71,059; 

but  in  no  case  less  than  zero. 

3.  Seaside  Area 

The  data  available  and  the  amount  of  water  drawn  from 
the  Seaside  wells  did  not  justify  a similar  model  of  the 
Seaside  aquifer.  The  Zone  11  investigation  places  safe  yield 
figure  for  Seaside  at  2,200  acre  feet  per  year  [Ref.  2]. 

This  study  will  treat  that  as  a constant  annual  import  of 
water  to  the  peninsula. 


IV.  THE  OBJECTIVE  FUNCTION  - AN  EXPERIMENTAL 
DETERMINATION  OF  UTILITY 


The  objective  of  this  model  is  directly  related  to  the 
value  that  water  has  to  the  residents.  A naive  approach 
would  be  to  assume  that  twice  as  much  water  should  provide 
twice  as  much  value.  If  this  were  the  case  it  would  not  be 
hard  to  see  that  the  solution  to  the  problem  is  to  consume 
all  water  as  soon  as  it  becomes  available.  But  is  this  a 
reasonable  argument?  It  is  equivalent  to  saying  that  out  of 
100  gallons  of  water  made  available,  the  second  fifty  gallons 
provides  the  same  amount  of  value  that  would  have  been  pro- 
vided by  only  having  50  gallons  to  begin  with.  Most  will 
agree  that  this  is  not  true;  that  the  value  of  water  per  unit 
volume  decreases  as  total  amount  of  water  available  increases 
Hence  the  need  for  a utility  function  to  provide  an  index  of 
how  the  value  of  water  changes  as  more  becomes  available. 
Reference  15  presents  an  excellent  review  of  utility  theory 
and  experimental  procedures  for  obtaining  information  about 
individual  utility.  This  study  assumes  that  a utility  func- 
tion for  water  exists  for  residents  of  the  Monterey  Peninsula 
What  the  model  needs  is  an  indexing  function  that  will  assign 
to  the  amount  of  water  D some  index  U(D)  that  will  permit 
comparisons  in  utilities  of  various  alternatives. 

Since  the  origin  and  unit  of  the  utility  scale  can  be 
assigned  arbitrarily,  one  can  select  any  2 alternatives  from 
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a given  set  of  alternatives  and  assign  fixed  values  of  utility 
to  them.  Then,  through  the  use  of  proposed  "lotteries" 
involving  those  2 alternatives,  individuals  may  be  induced 
to  reveal  on  the  same  utility  scale,  their  utility  of  some 
of  the  other  alternatives  available. 

An  experiment  based  upon  this  concept  was  conducted  in 
the  form  of  a survey  among  several  faculty  members  of  the 
Naval  Postgraduate  School  who  were  also  residents  of  the 
Monterey  Peninsula. 

A set  of  alternatives  was  presented  to  each  member  ques- 
tioned. Each  alternative  represented  a fixed  per  person 
daily  consumption  limit  to  which  the  population  would  be 
restricted  for  an  extended  period.  The  alternatives  ranged 
on  a continuous  scale  from  a strict  rationing  value  of  50 
gallons  per  person  per  day  (gpd)  to  unrestricted  consumption 
(arbitrarily  set  at  200  gpd) . 

Four  different  lotteries  were  then  presented  using  mix- 
tures of  the  2 extreme  alternatives  50  and  200.  Those 
questioned  were  asked  to  provide  the  single  alternative  that 
they  would  consider  equivalent  to  the  given  lottery. 

Computation  of  utilities  was  accomplished  by  fixing  the 
utility  of  50  gpd  at  0 and  the  utility  of  200  gpd  at  100.  For 
each  lottery  such  that  50  gpd  occurred  with  probability  p and 
200  with  probability  1-p,  the  utility  of  the  lottery  is  simply 
pU(50)  + (1-p)  U(200)  = (1-p)  (100)  . Thus  a utility  value  for 
each  answer  given  could  be  assigned  on  the  basis  of  the 
equivalent  lottery. 
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The  particular  lotteries  presented  yielded  those  alterna- 
tives whose  utilities  were  50,  60,  80  and  90  which  together 
with  0 and  100  already  fixed  gave  6 data  points  with  which 
to  estimate  the  utility  function.  The  resulting  utility 
function  was  constructed  using  the  median  responses  to  the 
survey  to  reduce  the  influence  of  poorly  answered 
questionaires . 

Figure  7 shows  the  points  that  resulted  from  the  survey. 
The  curve  in  the  figure  represents  a close  smooth  approxima- 


tion to  the  utility  function.  The  equation  of  the  curve  is 
9 


given  by: 


r 


U(D) 


(1. 617) D - 80.105 


for  D < 65  gpd 


(79.154)  /D  - (37.129)  Ln  D for  65  <_  D <_  175  gpd 


- ( 2 . 72 ) D - 281.354 


( . 0595) D + 87.83 


for  175  < D < 200  gpd 


A complete  copy  of  the  survey  is  provided  in  Appendix  B. 


a Empirical  data  was  created  from  a smooth  fit  drawn  through 
the  survey  points  and  multiple  regression  applied  to  the  sum 
of  appropriately  shaped  curves  (parabola,  line  and  logarithm) . 
Both  ends  were  then  refined  by  attaching  smooth  linear  extensions 


38 


V.  EXERCISING  THE  MODEL 


A.  THE  MODEL 

Coding  of  the  model  was  done  in  Fortran  IV  programming 
language  for  use  on  an  IBM  360/67  computer.  The  output  was 
designed  to  provide  a chart  which  the  decision  maker  could 
enter  with  a 1 May  state  of  the  system  and  extract  the  opti- 
mal decision  for  the  next  year's  consumption.  A copy  of  the 
program  can  be  found  at  the  end  of  this  study. 

The  model  itself  is  discrete  with  respect  to  the  storage 
supply,  the  feasible  consumption  decisions  and  rainfall. 
Storage  values  are  described  in  1,000  acre  feet  increments 
up  to  the  maximum  capacity,  each  with  three  possible  condi- 
tions of  wetness.  Thus  for  a maximum  capacity  of  38,000  acre 
feet,  there  are  117  possible  states  for  entry  into  the  chart. 
The  permissible  decisions  for  each  state  are  also  in  1,000 
acre  feet  increments  ranging  from  0 to  the  amount  in  storage 
indicated  by  that  state.  Rainfall  was  likewise  made  discrete 
to  facilitate  computation  of  expected  values  of  utility. 

Utilities  are  computed  using  the  function  derived  earlier 
as  typical  for  each  resident  of  the  peninsula.  Since  the 
utility  function  is  dependent  upon  net  consumption  by  an  indi- 
vidual, a conversion  is  required  from  the  decision  in  acre 
feet  to  the  number  of  gallons  per  person  per  day  that  it 
represents.  Approximately  5%  of  all  production  is  lost  to 
waste  before  metering  by  Cal  Am  and  only  86%  of  metered 
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production  goes  to  residential  use.  In  addition  the  model 
assumes  an  automatic  importation  of  2,200  acre  feet  from  the 
Seaside  aquifer.  The  current  population  figure  used  by  the 
model  is  91,000  residents  serviced  by  Cal  Am.  Therefore  for 
a given  decision  D of  acre  feet  to  consume,  only 
(365)  (9l'000)  (-95M-86>  (D-2840+2200)  gpd  is  available  for 

consumption;  where  7.5  = # gal  per  cubic  ft,  43,560  = # 
cubic  feet  per  acre  feet  and  2840  = # acre  ft  of  private  well 
extraction  and  subsurface  losses  which  are  assumed  uncontrollable. 

B.  ASSUMPTIONS 

The  following  assumptions  were  made  to  facilitate  the 
running  of  the  model.  Some  are  obviously  not  true  in  an  exact 
sense  but  at  the  same  time  they  are  not  expected  to  disrupt 
the  accuracy  to  any  significant  degree. 

Assumption  1:  Runoff  is  assumed  to  recharge  the  aquifer 
instantly.  That  is,  if  there  is  sufficient  capacity  left  in 
the  alluvium  to  hold  the  computed  runoff,  then  it  will  all 
percolate  into  the  aquifer  before  reaching  the  ocean. 

Assumption  2:  Given  the  opportunity  to  consume  water,  the 
population  will  consume  it.  This  assumption  is  made  more 
realistic  by  bounding  production  by  the  maximum  historical 
figures  from  recent  normal  consumption.  The  model  essentially 
stops  considering  any  decisions  that  would  provide  more  extrac- 
tion than  the  population  would  freely  consume.  Figure  8 shows 
a normal  consumption  pattern  under  no  restrictions.  A produc- 
tion limit  of  20,000  acre  feet  is  placed  on  the  model  to 


preclude  average  daily  per  capita  consumption  that  would 
exceed  that  pattern. 

Assumption  3:  It  is  assumed  that  the  California  American 
Water  Company  has  the  capability  of  meeting  the  production 
demands  recommended  by  the  model.  That  is,  given  that  .the 
model  recommends  a certain  figure  for  optimal  consumption, 
the  water  company  will  be  able  to  produce  that  quantity. 

The  1975-1977  drought  revealed  that  this  is  in  fact  not  the 
case.  This  problem  will  be  discussed  in  the  last  section. 

C . RESULTS 

The  output  of  the  model  is  presented  following  Appendix  B. 
All  possible  states  are  represented.  There  are  3 pages  of 
output,  each  representing  a separate  initial  wetness  condi- 
tion. One  selects  dry,  normal  or  wet  as  a result  of  the 
current  estimate  of  annual  runoff  due  to  current  rainfall 
(since  stream  flow  may  not  be  complete,  this  can  be  estimated 
using  the  runoff  model  of  Section  III) . Find  the  column 
whose  first  row  entry  equals  the  current  estimated  state  and 
read  the  optimal  decision  in  acre  feet  in  column  2 and  the 
optimal  rationing  policy  in  column  4.  Column  3 contains  the 
value  of  the  objective  function  due  to  the  particular  optimal 


decision. 


VI.  CONCLUSIONS  AND  FURTHER  OBSERVATIONS 


A.  CONCLUSIONS 

A careful  study  of  the  results  of  the  model  reveal  that 
the  optimal  decision  is  always  to  consume  all  available  supply 
when  that  supply  falls  below  the  normal  unrestricted  consump- 
tion established  by  past  experience.  The  result  is  not  sur- 
prising in  view  of  the  normal  hydrologic  cycle  for  the  area. 

With  a virtual  capacity  conservatively  estimated  at  38,000 
acre  feet  and  average  annual  runoff  in  excess  of  60,000  acre 
feet,  the  model  has  reason  to  be  optimistic  about  future 
replenishment  of  supplies.  This  optimism  coupled  with  the 
rather  mild  rate  of  change  in  the  slope  of  the  utility  func- 
tion is  responsible  for  the  conclusion  reached. 

This  does  not  imply  that  there  is  no  need  for  rationing 
during  consecutive  years  of  drought.  It  merely  implies  that 
under  the  conditions  of  the  model,  there  is  nothing  to  be 
gained  by  sacrificing  present  consumption  to  add  to  future 
supplies.  If  current  supplies  are  low,  of  course,  rationing 
is  called  for  to  preclude  running  out  in  the  current  year. 

B . FURTHER  OBSERVATIONS 

An  interesting  fact  was  brought  out  by  the  model.  According 
to  records,  the  1974-75  runoff  was  above  average.  Assuming 
38,000  acre  feet  in  storage  on  1 May  1975  and  computing  forward 
from  the  data  available,  the  amount  of  storage  on  1 May,  1977 


was  estimated  to  have  been  about  16,500  acre  feet.  The  model 
recommends  consuming  all  which  equates  to  a rate  of  125  gpd 
for  the  year.  Why  then  did  Cal  Am  impose  strict  rationing  of 
50  gpd? 

The  answer  to  this  question  requires  re-examination  of 
the  extent  of  the  Carmel  Valley  alluvium.  Figure  9 shows  a 
vertical  cross-section  of  the  aquifer  along  with  the  boun- 
daries of  the  Cal  Am  well  fields  and  elevation  of  the  water 
table  during  and  after  the  drought. 

Since  supplies  were  considerably  below  normal  in  early 
winter  of  1977,  the  major  source  of  production  by  Cal  Am 
would  necessarily  come  from  their  wells.  Normal  well  produc- 
tion was  fully  restored  in  Feb.  1978  at  12,760  gal  per  minute. 
This  compares  with  the  early  December  1977  maximum  rate  of 
1630  gal  per  minute.  Figure  9 illustrates  the  reason  for  the 
drastic  reduction  in  extraction  capability. 

The  December  1977  level  of  the  water  table,  although  not 
excessively  low  throughout,  was  less  than  50%  of  normal  levels 
in  the  region  of  the  aquifer  containing  the  Cal  Am  well  fields. 
As  a result,  production  capacity  of  the  Carmel  Valley  alluvium 
was  limited  to  less  than  3,000  acre  feet  per  year.  Anticipa- 
tion of  only  3,000  acre  feet  plus  best  estimates  of  an  addi- 
tional 4,000-5,000  acre  feet  of  combined  production  from 
Seaside  and  the  reservoirs  prompted  the  decision  by  Cal  Am. 

The  obvious  solution  to  this  problem  after  observing 
Figure  9 is  to  place  wells  in  the  western  most  portion  of  the 
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CROSS-SECTION  OF  C.V.  ALLUVIUM  (NOT  TO  SCALE) 


FIGURE  9 


I 


aquifer.  This  has  been  proposed  and  will  most  likely  be 
resolved  in  the  near  future.10 

A final  observation  is  offered  concerning  the  controver- 
sial topic  of  the  future  growth  of  the  Carmel  Valley/Monterey 
peninsula  area;  in  particular  how  the  increased  demands  for 
water  will  be  met.  -According  to  the  runoff  model  derived 
earlier,  water  needs  of  the  future  will  be  supplied  in  great 
quantity.  The  deficiency  lies  in  adequate  storage.  At 
present  an  average  of  60,000  acre  feet  of  stream  flow  is 
lost  to  sea  each  year,  runoff  which  if  stored  could  provide 
more  than  twice  the  present  demand  for  water. 


There  are  questions  of  water  quality  that  enter  here 
which  although  solvable  do  have  an  impact  on  questions  of 
where  to  locate  new  wells. 
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APPENDIX  A 


MATHEMATICAL  REPRESENTATION  OF  THE  MODEL 


Let: 


n = subscript  indicating  the  number  of  years 
remaining  in  the  time  horizon  selected 
for  optimization 

a « discount  factor  0 < a < 1 

X = (x  ,w  ) = State  vector  composed  of 
n n n 

xn  = amount  of  water  in  storage  and 

w = wetness  condition  of  the  soil 
n 

r = total  annual  rainfall,  a random  variable 
n 

Vn(Xn)  = Optimal  value  of  the  objective  function 
given  the  state  X 

n 

U(D)  = utility  of  Dgpd  of  water  consumption 

D^  = a feasible  consumption  decision 

C = maximum  storage  capacity  for  water 

R(r,w)  = runoff  due  to  r inches  of  rainfall  under 
conditions  of  wetness  = w 


The  model  was  solved  recursively  using  dynamic  programming 
techniques  as  follows: 

Let 


V (X)  = 0 (there  is  no  utility  if  there  is 

° no  time) 
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bow , 


then 


V (X,)  * max  (U (D . ) + aE[V  (X_)J  - max  [U(D.)] 

-► 

and  given  V (X  ) , the  model  computes  V (X  ^.)  as  follows: 
n n n+i  n+i 


Vn+1 IXn+l> 


°-Di— xn+l 


{U (D; ) + aE  (V  (X  ) ] } 
i n n 


where 


X = min[C;X  , - D.  + R(r  ,w  )] 
n n+l  i n n 


E[S]  "Expected  Value"  of  S. 


The  values  of  parameters  used  were: 


C = 38,000  acre  feet 


R(r,w)  = max{5296r  + 8433w  - 71,058;0) 


1 if  R(r  . ,w  ,)  < 1160  acre  feet 

n-1  n- 1 

2 if  1160  < R(r  . ,w  .)  < 148793  acre  feet 

— n-l  n-i  — 

3 if  148,793  < R(r  . ,wb  . ) 

n-l  n-i 
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The  selection  of  a value  for  n can  be  somewhat  arbitrary. 
Since  future  utilities  are  discounted,  there  should  be  some 
value  of  n beyond  which  the  solution  becomes  independent  of 


50 


APPENDIX  B 


THE  SURVEY 


The  following  survey  is  being  conducted  to  determine  the 
value  of  water  to  the  population  of  the  Monterey  Peninsula. 
More  specifically,  as  the  amount  of  water  for  consumption 
increases , what  is  the  relative  increase  in  value  to  the  con- 
sumer? The  results  will  be  used  in  a thesis  study  at  the 
operations  research  dept,  of  NPS. 

The  basic  unit  of  measure  for  water  consumption  will  be 
expressed  in  gallons  per  person  per  day  (gpd) ; however, 
compliance  with  consumption  restrictions  will  be  monitored 
monthly.  Consumption  restrictions  (rationing)  are  set  each 
year  for  the  entire  year  and  can  range  from  a minimum  (strict) 
of  50  gpd  on  up  to  200  gpd  (considered  unrestricted) . 

Consider  the  time  period  from  now  through  the  next  10 
years.  You  will  be  presented  with  an  expected  rationing 
pattern  over  those  ten  years  that  evenly  distributes  a given 
number  of  strict  rationing  years  among  years  with  no  consump- 
tion restrictions  whatsoever.  On  the  right  please  enter  the 
fixed  rationing  policy  that  you  consider  to  be  equivalent  in 
value  for  you. 


For  Example : 


Policy  A 

8 years  at  50  gpd  and 
2 years  unrestricted 


or 


Policy  B 

10  years  of  continuous 
rationing  fixed  at 
60  gpd. 


Indicates  that  you  are  indifferent  between  rationing  at  60  gpd 
100  percent  of  the  time  and  rationing  at  50  gpd  80  percent  of 
the  time. 

Attached  to  the  survey  you  will  find  some  relevant 
statistics  and  current  usage  information  that  should  be  of 
considerable  value  in  determining  your  choices. 

Upon  completion  please  fold  and  place  the  survey  in  the 
guard  mail  for  forwarding.  Thank  you. 
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QUESTIONNAIRE 


I.  For  this  section  please  enter  the  rationing  value  in 
Policy  B which  will  cause  you  to  be  indifferent  between 
Policy  A and  Policy  B. 


Policy  A 

Policy  B 

1. 

5-years  of  strict  rationing 
at  50  gpd  and  5 years  of 
unrestricted  consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 

Policy  A 

Policy  B 

2. 

4 years  of  strict  rationing 
at  50  gpd  and  6 years 
unrestricted  consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 

Policy  A 

Policy  B 

3. 

2 years  of  strict  rationing 
at  50  gpd  and  8 years 
unrestricted  consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 

Policy  A 

Policy  B 

4. 

1 year  of  strict  rationing 
at  50  gpd  and  9 years 
unrestricted  'consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 

II.  For  this  section  please  enter 
majority  of  local  residents  would 

the  values  that  you  feel  the 
enter  if  asked. 

Policy  A 

Policy  B 

1. 

5-years  of  strict  rationing 
at  50  gpd  and  5 years  of 
unrestricted  consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 

Policy  A 

Policy  B 

2. 

4 years  of  strict  rationing 
at  50  gpd  and  6 years 
unrestricted  consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 

Policy  A 

Policy  B 

3. 

2 years  of  strict  rationing 
at  50  gpd  and  8 years 
unrestricted  consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 

Policy  A 

Policy  B 

4. 

1 year  of  strict  rationing 
at  50  gpd  and  9 years 
unrestricted  consumption 

OR 

10  years  of  continuous 
rationing  fixed  at 
gpd 
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I.  STATISTICS 

CONSUMPTIVE  ACTION  QUANTITY  OR  RATE  OF  CONSUMPTION 


A. 

FILLING  AVERAGE  TUB 

30 

Gals 

B. 

FLUSHING  TOILET 

5 

Gals 

C. 

5 MINUTE  SHOWER 

30  Gals  (6  Gals/mi: 

D. 

DISHWASHER  (ONE  LOAD) 

15 

Gals 

E 

AUTOMATIC  CLOTHES  WASHER 

1 . MINIMUM  LOAD 

25 

Gals 

2 . MAXIMUM  LOAD 

50 

Gf  Is 

F. 

KITCHEN  FAUCET 

6 

Gals/min 

G. 

COMMON  GARDEN  HOSE 

250 

Gals/hr 

H. 

COMMON  LAWN  SPRINKLER 

120 

Gals/hr 

II.  Personal  Hygiene  requirements  for  an  adult  range  from  an 
estimated  minimum  of  21  gals/day  to  an  extravagant  maximum  in 
excess  of  75  gals/day.  A reasonable  figure  under  normal 
conditions  is  about  45  gals/day.  This  figure  accounts  for 
drinking,  bathing  and  hygiene  only. 

III.  Consumption  data  from  Cal-Am  Water  Company  indicates  that 
during  the  normal  years  preceeding  the  1975-77  drought,  the 
consumption  of  water  by  local  residents  varied  from  a rainy 
winter  season  low  rate  of  about  85  gpd  to  a late  summer  high 
of  about  175  gpd.  The  difference  is  due  primarily  to  outdoor 
uses  of  water.  Particularly  maintenance  of  lawns  and  gardens. 

For  example,  the  recommended  rate  of  water  application  for 
a healthy  lawn  is  1 inch  per  week  for  the  coastal  slopes  and 
1.5  inches  per  week  for  the  inland  coastal  areas  of  central 
California. 

Note ; 1 in/wk  on  5000  sq  ft  (~l/8  acre)  = 3125  gal/wk  = 446 
gal/day . This  equates  to  over  110  gpd  for  a family  of  four. 

To  prevent  a lawn  from  dying  during  drought  conditions, 
one  is  recommended  to  water  for  a minimum  of  15  minutes  per 
week.  For  the  same  5000  sq  ft  of  lawn,  this  amounts  to  105 
gals/week  or  15  gals/day. 

Recall  that  gpd  refers  to  gals  per  person  per  day. 
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RESULTS  WHEN  INITIAL  CONDITIONS  ARE  WET 
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PROGRAM  LISTING 


COMMON  VN(3,3,  70  ) , PROOUC.  C AP  AC,  POPUL 
EXTERNAL  UTILS,  RUNOFF,  CONVER 
INTEGER  01  MEN 

fOPUL -91000. 

AP AC-38000. 

PROOUC -20000. 

OIMEN-IFIX(CAP AC/1000. *1. 5) 

CALL  EXECUT(DIMEN) 

STOP  ' 

END 

SUBROUTINE  EXECUT(OIMEN) 

COMMON  VN( 3,  3,  70 ), PROOUC,  C APAC,  POPUL 
INTEGER  DIMEN 

DIMENSION  STORG!  70 ) , STATE ( 3 ) , GPD( 70 ) 

REAL *8  STATE  /' DRY ', 'NORMAL ', 'WET ' / 

L-0 

DO  SO  1-1,  OIMEN 

STORG(  I ) -FLOAT!  1000-0-1)) 

SO  CONTINUE 

00  1 1-1,  3 

DO  2 K-l,  OIMEN 

VN(  I , 1,  K)  ■ AM  I Nl( STORG! K), PROOUC) 

VN ( I, 2, K)  -UTILS(VN( I, 1,  K) ) 

VN( I , 3, K) -VN( I , 2, X) 

2 CONTINUE 

1 CONTINUE 

7 ICEY-0 

. DO  3 1-1, 3 

CALL  NXTSTG! I,  OIMEN) 

3 CONTINUE 

L-L-l 

00  4 I -1,  3 

DO  S X -1,  OIMEN 

IF((VN(I,2,K)-VN(I,3,  K)).GT.10.  ) ICEY-1 
VN(  I,  3,  X)  -VN!  I,  2,10 
S CONTINUE 

4 CONTINUE 

IF(L. EQ.10)  GO  TO  20 
IFOCEY.  EQ.l)  GO  TO  7 
20  DO  9 INOEX-1,  3 

00  10  I I -1, OIMEN 

GPO( I I ) -CONVER! VN( INDEX. 1,  II)) 

10  CONTINUE 

WRITE(6, 100)  STATE! I NOEX) 

100  FORMAT!////'  RESULTS  WHEN  INITIAL  CONDITIONS  ARE  ',  A6) 

DO  8 Ml-1,  OIMEN,  12 
M2-M1-11 

WRITE (6, 101)  (STORG! I ), I -Ml, M2), (VN! INDEX,  1,  J),  J-M1.M2), !VN! 
1INDEX,2,K),K4U,M2),  (GPO(L),  l -Ml,  M2) 

101  FORMAT!//'  WATER  ON  HANO:  ' , 12 ( IX,  F7.  0 ) / ' OPTIMAL  CONSUMP:  ',1 

12! IX, F7. 0)/ ' TOTAL  UTILITY:  ', 12! IX, F7. 0)/ 

2'  GAL/PERSON/DAY:  ' , 12! IX, F7. 0 ) // ) 

8 CONTINUE 

9 CONTINUE 

RETURN 
ENO 
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SUBROUTINE  NXTSTG( I WET, DI  MEN) 

COMMON  VN( 3, 3,  70 ) , PRODUC, CAP AC, POPUL 
INTEGER  DIMEN 
01  MENS  I ON  RAINFU36) 

DATA  RAINFl  /.  002, 3*.  0005,  4* . 0006, . 009,  2*. Oil,  . 015, . 016, . 024, . 035, 
1.  047,  . 078,.  10  2,.  12,.  06,.  04,.  035,  . 025,  2*.  02,  5*.  015,  2*.  01 2,.  Oil, 

22*. 01, .169/ 

ALPHA-.  9 
DO  1 1 1-1,  D I MEN 
OMAX-O. 

UMAX— 10000. 

STORGl-FLOATC  1000* (I  1-1) ) 

PUMPGE -AMI N1 ( STORG1,  PRODUC) 

M AXD  - I F I X(  PUMPGE/  1000. *1.  5) 

DO  2 12-1, MAXD 

DEC  I S-FLOAT(  1000*(  I 2-1) ) 

SUM-0. 

00  3 13-1,  36 
RAI N-FLO AT ( I 3) 

RETURN-RUNOFF ( RAI N.  IWET) 

STORG2- AM  INK  STORGl  ♦RETURN-DEC  I S,  CAP  AC) 

1 WET 2 -2 

IFtRETURN. GT. 148793. ) IWET2-3 
IF(RETURN.  LT. 11601.  ) IWET2-1 
INDEX -I  FI X(STORG2/ 1000.  *1.  5) 

SUM -SUM-RAI NFL ( I 3)*VN( IWET2,  3,  INDEX) 

3 CONTINUE 

UEXP-  UTILS(DECIS) ♦ALPHA* SUM 
I F (UEXP.  LT. UM AX)  GO  TO  2 
UMAX -UEXP 
DMAX -DEC  IS 
2 CONTINUE 

VN< IWET, 1,  I 1) -OMAX 
VN(  IWET,  2,  ID -UMAX 
1 CONTINUE 

RETURN 
END 
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FUNCTION  RUNOFFCRAIN, IWET) 

RUNOFF >5296. 38932*RAIN*8  433.  159*FLOAT< I WET) -71058.  62 
RUNOFF-AMAX1(RUNOFF,0. ) 

RETURN 

END 


FUNCTION  UTILS(WATER) 

GPD*CONVER( WATER) 

IF(GPO. LT.65. ) UTILS-1. 617«GPD-80.  105 

I F(GPO. GE. 65 . . AND.GPO. LE.175.  ) UTI LS-(  79.  1S4)*SQRT(GPD) -<37.  129  )* 
lALOGlGPD)-2.  72*GPD-231.  354 
I F(GPO. GT. 175 ) UTIlS-(. Q59S7)*GPD*37. 83 
RETURN 
END 


FUNCTION  CONVER(X) 

COMMON  VN( 3,  3,  70 ) , PRODUC, CAPAC,  POPUL 
CONVER*( 731. 27 ) * ( X-640. )/POPUL 
CONVER-AMAXHCONVER,  0.  ) 

RETURN 

END 
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